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ABSTRACI 


A  ccuputer  model  was  developed  to  simulate  a 
projectile  ram- propellant  feed  system  for  a  Liquid 
Propellant  Gun.  Using  a  lumped  parameter  approach,  a 
set  of  sinultaneous  differential  equations  was  derived 
for  the  complex  interaction  of  the  propellant  fluid, 
the  driving  injector  and  the  projectile.  The  computer 
model  was  verified  against  a  20  mm  experimental 
apparatus.  Injector  displacement,  projectile 
displacement,  and  chamber  pressure  were  compared  for  a 
nominal  driving  pressure  of  140  psi.  The  important 
system  parameters  affecting  projectile  ram  time  and 
chamber  pressure  cscillations  were  investigated  and 
potential  problem  areas  for  testing  with  actual 
propellant  vere  identified. 
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!•  INTRODUCTION 


As  the  state  cf  the  art  of  weapons  technology  improves, 
mere  complex  systems  are  continually  being  developed.  To 
meet  future  needs  these  systems  must  have  improved 
performance,  less  weight  and  volume,  less  complicated 
logistics,  and  at  the  same  time  be  more  cost  effective  than 
current  systems.  One  candidate  for  these  future  veapon 
systems  is  the  liquid  propellant  gun  (LPG) .  Numerous 
invest igatiers  over  the  past  thirty  years  have  demonstrated 
its  potential  if  the  technology  can  be  developed  to  utilize 
a  liquid  propellant  in  a  large  caliber  gun. 

One  important  area  that  needs  investigation  is  the 
aspect  of  fluid  prcpellant  handling.  The  ordnance  designer 
must  kncii  the  necessary  propellant  flow  rates  and  pressure 
regimes  needed  to  achieve  the  desired  system  performance. 
These  parameters  must  be  within  allowable  safety  margins  for 
the  handling  and  use  of  the  explosive  propellant.  One 
possible  design  of  an  LPG  which  would  allow  a  very  high 
firing  rate  utilizes  the  propellant  to  ram  the  loaded 
projectile.  Combining  the  normally  separate  propellant  load 
and  ram  cycles  decreases  the  required  movement  of  mechanical 
parts.  Thus,  using  this  method,  it  is  feasible  that  firing 
rates  in  excess  cf  fpur  to  five  times  present  rates  could  be 
achieved.  This  thesis  was  directed  toward  identifying  the 
important  fluid  dynamic  parameters  involved  with  such  a 
projectile  ram  feed  system. 
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II.  BACKGBOONC 


During  1S74  and  early  1975,  an  investigation  was 
conducted  at  the  Naval  Postgraduate  School  to  study,  both 
analytically  and  experimentally,  the  fluid  dynamics  of  a 
liguid  propellant  under  conditions  similar  to  those  which 
would  exist  in  a  rapid-fire  LPG  feed  system.  The  results  of 
this  investigation  were  to  be  used  to  establish  such  LPG 
design  and  performance  parameters  as  time-to-load,  injection 
supply  pressure,  injection  system  configuration,  ullage, 
charge-tc-mass  ratio,  caliber  size,  and  projectile  mass.  It 
was  hoped  that  this  investigation  would  identify  any 
potential  problem  areas  for  further  detailed  research. 

The  basic  objective  of  the  experimental  portion  cf  the 
investigaticii  was  to  identify  what  fluid  dynamic 
characteristics  of  a  liquid  propellant  feed  system  would 
limit  leading  times  and  hence  rates  of  fire.  A  20  mm 
experimental  model  of  a  basic  propellant  feed  system  was 
designed  and  built.'  Data  on  injector  displacement,  breech 
chamber  pressure,  and  ram  gas  pressure  (input  driving 
pressure) ,  were  recorded  for  driving  pressures  between  50 
and  220  psig  in  10  psi  increments.  As  reported  in  Ref.  1, 
it  was  found  that  the  instantaneous  behavior  of  the  chamber 
pressure  was  the  result  of  a  complex  interaction  of  inertia 
forces,  visccus  forces,  and  the  unsteady  motion  of  the 
fluid.  The  experiments  demonstrated  that  frictional  and 
inertial  effects  were  significant  during  the  movements  of 
the  injector  and  the  projectile  slug.  Once  the  projectile 
slug  stopped,  the  effect  of  entrapped  gas  in  the  fluid 
caused  large  breech  chamber  pressure  oscillations.  On 
several  runs,  sub-atmospheric  pressures  were  experienced 
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which  suggested  the  possibility  of  cavitation  and  hence 
vapor-phase  igniticn-  Over  all^  it  was  found  that  the  ram 
tine  had  a  guadratic  dependence  upon  ram  pressure. 

The  analytical  portion  of  the  investigation  was  directed 
toward  predicting  the  pressure  and  flew  rate  of  the 
propellant  atd  the  projectile  slug  motion  during  an  LPG 
loading  cycle.  The  mathematical  model  which  was  developed 
and  programmed  on  an  analog  computer  predicted  the  position, 
velocity,  and  acceleration  of  the  projectile  and  the 
pressure  at  various  points  in  the  system  as  functions  of 
time.  The  model  was  tested  against  experimental  results  and 
found  to  he  adequate  for  the  prediciton  of  projectile  ram 
time.  An  analysis  was  performed  using  the  model  to  indicate 
the  areas  of  system  redesign  likely  to  be  most  profitable 
and  to  obtain  preliminary  prediotions  of  IFG  loading  system 
performance  under  a  variety  of  design  conditions.  The 
results  of  the  analytical  study,  as  well  as  a  summary  of  the 
experimental  work  is  contained  in  Bef.  2. 

The  aforementioned  model  was  the  starting  point  for  the 
present  studj.  To  better  understand  the  derivations  in  the 
following  sections.  Figure  (1)  indicates  the  geometry  of  a 
hasic  projectile  ram  propellant  feed  system.  An  injector 
chamber  is  filled  with  propellant,  which  is  then  pumped  into 
the  gun  hreech  by  applying  ram  pressure  to  the  ram  side  of 
the  injector  piston.  This  can  be  accomplished  by  using  high 
pressure  gas  from  an  accumulator  or  by  using  an  hydraulic 
drive  sjstei.  The  force  exerted  by  the  injector  piston 
drives  the  fluid  propellant  through  the  connecting  line  into 
the  gun  breech.  The  rapidly  accelerated  propellant  drives 
the  already  loaded  projectile  to  its  seated  position  ready 
for  firing. 
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Figure  1  -  BASIC  GEOMETBI 
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Ihe  cccflex  nature  of  this  dynamic  ram  process  is  not 
easy  to  describe  with  mathematical  eguations.  Hany 
researchers  have  approached  the  problem  using  wave  mechanics 
and  partial  differential  equations.  These  equations  were 
then  modified  according  to  the  propagation  characteristics 
of  the  system  and  the  various  boundary  conditions 
encountered  therein.  This  usually  led  to  involved  finite 
difference  methods  of  solution  on  a  digital  computer. 

Another  possible  approach  considers  the  kinetic  energy 
involved  in  the  ram  process.  This  is  the  approach  that  was 
used  to  derive  the  governing  eguations  for  the  previously 
described  analog  computer  model.  Beginning  with  the  input 
side^  a  force  balance  on  the  injector  piston  yielded  an 
equation  fer  the  pressure  in  the  injector  chamber  in  terms 
of  the  injector  motion.  Writing  Bernoulli's  equations  for 
head  losses  at  the  ejepansion  to  the  breech  chamber ,  the 
ccntracticn  to  the  connecting  line,  and  an  orifice  in  the 
connecting  line  resulted  in  an  equation  for  the  pressure 
drop  in  the  connecting  line  in  terms  of  the  square  of  the 
fluid  velocity  and  tie  fluid  acceleration.  A  force  balance 
on  the  projectile  slug  yielded  a  third  equation  which 
described  the  breech  chamber  pressure  in  terms  of  the  slug 
motion.  To  solve  these  equations,  some  method  of  relating 
the  motions  in  the  injector,  connecting  line,  and  breech 
chamber  had  to  be  established.  It  was  assumed  that  the 
fluid  was  incompressible  and  therefore  these  motions  were 
identical.  Tc  facilitate  the  analog  simulation  it  was  also 
assumed  that  the  irput  ram  pressure  was  a  step  input.  In 
the  process  of  converting  these  equations  to  an  integral 
form  for  wiring  on  an  analog  computer,  all  variables  were 
normalized  tc  permit  scaling. 

Despite  these  limiting  assumptions,  the  analog  computer 
model  was  able  to  predict  the  time  to  ram  the  projectile  and 
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the  chamter  pressures  which  cccur  duriog  the  time  that  the 
injector  is  in  motion.  Due  to  the  incompressitility 
assumption,  the  model  was  unable  to  account  for  the  pressure 
transients  or  "water  hammer"  type  pressure  oscillations 
which  occurred  when  the  projectile  was  seated  and  the 
propellant  fluid  was  being  decelerated.  These  pressure 
oscillations  are  important  to  the  system  designer  because 
the  peak  breech  chamber  load  pressure  is  experienced  during 
these  oscillations.  It  is  also  possible  that  these 
oscillations  could  interfere  with  the  uniformity  cf  the 
ignition  and  subseguent  combustion  of  the  propellant. 

£  straight  fcmard  attempt  was  made  to  extend  this 
analog  mcdel  by  first  converting  it  to  a  digital  conputer 
pregram.  This  allpwed  an  accurate  representation  cf  the 
input  ram  pressure  and  the  use  of  unnormalized  variables  to 
be  incorporated.  These  improvements  increased  the  accuracy 
of  the  mcdel  but  the  effect  of  the  compressibility  of  the 
propellant  fluid  still  was  not  taken  into  account.  Trying 
to  add  compressibility  effects  to  this  model  by  adding  time 
derivatives  cf  the  pressure  terms  to  the  governing  equations 
unlinked  the  motions  of  the  injector,  the  connecting  line 
fluid,  and  the  projectile  slug  causing  a  problem  with  too 
many  unknown  variables  for  the  number  of  equations  involved. 

It  was  decided  that  a  new  approach  should  be  tried  which 
would  adequately  describe  the  system  pressure  oscillations. 
The  approach  would  feature  an  engineering  mcdel  which  would 
be  easy  to  use,  adaptable  to  any  LPG  feed  system,  and  not 
obscure  the  interaction  of  system  parameters  by  complicated 
mathematics. 
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III.  DEBIVATION  Cf  MODEL 


A.  LDMPEI  tiEAMETEB  APPROACH 

The  apprcach  that  showed  the  greatest  promise  for 
modeling  the  LPG  feed  system  was  that  of  the  fluid 
transmission  line  cpncept.  This  approach^  which  has  tecome 
popular  in  the  last  five  to  ten  years,  is  based  on  a 
pressure-voltage  and  flow  velocity-current  analogy  with 
Electrical  Engineering  determinations.  It  is  an  outgrowth 
of  the  large  amount  of  effort  that  has  been  devoted  to 
investigating  fluid  line  transients.  Reference  3  is  a  good 
survey  of  this  field. 

For  complex  systems  this  approach  is  usually  simplified 
ty  using  the  approximation  of  lumped  parameters.  The 
effects  of  fluid  inertia,  capacitance,  and  resistance  are 
"lumped"  and  considered  to  act  only  in  discreet  areas  of  the 
system.  This  results  in  a  reduction  of  the  unknown 
parameters  due  to  the  assumed  lack  of  interaction  of  the 
different  fluid  effects.  This  approximation,  however, 
results  in  the  necessity  of  using  several  empirical 
constants  which  must  be  determined  by  fitting  experimental 
data.  The  ordinary  differential  equations  that  are  derived 
from  this  method  can  be  solved  either  by  Laplace 
transfer nation  or  ty  computer  integration. . 

To  consider  the  effect  of  fluid  inertia  it  is  assumed 
that  only  pressure  and  inertia  forces  are  present  and  that 
compressibility  effects  are  negligible  in  the  volume  under 
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coDsideraticr .  Then^ 


IV 


is  th€  propellant  fluid  density,  L 


where  ^  is  the  propellant  fluid  density,  L  is 

characteristic  length  and  V  =V  =V  because  the  flew  i 

1  2 

considered  inccmpressible  for  this  building  block  (Fig.2A). 


To  consider  the  effect  of  fluid  capacitance,  it  is 
assumed  that  only  compressibility  effects  are  important,  and 
that  inertia  and  resistance  effects  may  be  neglected  in  the 
vcluite  under  consideration.  Therefore, 

L  dP 

V  -  V  - - 

1  2  K  dT 

where  K  is  the  effective  bulk  modulus  cf  the  propellant 
fluid,  L  is  a  characteristic  length  and  P^=P^=P  (Fig.  2B) . 

Because  fluid  resistance  can  be  affected  by  sc  many 
different  parameters#  it  is  impossible  to  write  a  general 
equation  describing  the  pressure  drop  due  to  fluid 
resistance.  It  is  best  to  treat  it  empirically  using  an 
experimentally  derived  figure.  Thus,  in  the  volume  under 
consideraticr , 

P  -  P  =  Rv  V 

1  2 

where  Bv  is  a  function  of  fluid  velocity  and  (^^9* 

2C) .  Changes  in  cross  sectional  area  can  be  accounted  for 
by  using  appropriate  area  ratios. 

These  results  specify  three  building  blocks  which  can  be 
cciibined  in  any  sequence  to  model  the  dynamic 
characteristics  of  a  system.  The  complexity  of  the  model 
can  be  increased  to  any  degree  necessary  by  including  more 
and  more  cemhinatiens  of  these  three  basic  building  blocks. 
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A)  Inertia 
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B)  Capacitance 


C)  Resistance 
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E.  IPG  SYSIEM  GOVERNING  EQUATIONS 


Several  different  combinaticns  of  building  blocks  were 
tried  in  attempting  to  model  the  experimental  LPG  system. 
The  experimental  results  showed  that  inertial  and  resistance 
blocks  should  be  included  in  the  model.  The  fact  that 
pressure  transients  pccurred  in  the  breech  chamber  indicated 
that  a  capacitance  block  should  also  be  included.  To 
ttinimize  the  use  of  computer  time,  it  was  decided  to  start 
with  the  most  simple  model  thought  adequate  and  work  toward 
more  complex  models  as  the  quality  of  the  computer  results 
dictated. 

The  first  attenpt  considered  that  ccmpressibilt y  effects 
would  dcttinate  in  the  injector  chamber  as  the  accelerating 
piston  interacted  with  the  propellant  fluid.  The  resistance 
effects  cf  changes  in  the  feedline  cross  sectional  diameter 
and  associated  orifices  were  modeled  as  one  resistance  block 
at  the  beginning  of  the  connecting  line.  Inertial  effects 
were  thought  to  dominate  in  the  connecting  line.  This  would 
be  particularly  true  in  large  scale  models  where  the  length 
of  the  connecting  line  would  be  very  large.  It  was  thought 
that  compressibility  effects  should  again  dominate  in  the 
breech  chamber  as  the  accelerating  propellant  fluid  drove 
the  projectile  dcwn  the  breech  chamber.  Results  in  this 
case  were  governed  by  inertial  effects  because  the 
calculated  ccnnecting  line  velocity  was  too  great.  This  led 
to  a  decrease  in  the  chamber  pressure  as  the  projectile 
accelerated  until  a  large  negative  chamber  pressure  existed. 
In  addition  to  this  deficiency,  the  computed  system  pressure 
cscillaticns  were  not  sufficiently  damped. 

To  try  tc  adjust  the  connecting  line  velocity,  another 
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resistance  tlcck  was  added.  As  will  be  seen  in  the  section 
discussing  the  parameters  which  affect  the  pressure 
oscillaticns y  the  ^alue  of  this  resistance  coefficient  is 
picporticral  to  the  amount  of  damping  in  the  pressure 
oscillations.  Even  with  a  value  which  gave  the  proper 
damping,  the  breech  chamber  pressure  was  dominated  ty  the 
pressure  loss  due  to  the  accelerating  fluid  in  the 
connecting  line. 

To  reduce  the  dominance  of  the  fluid  inertia  effects,  it 
was  decided  tc  restructure  the  model.  To  model  the  injector 
chamber  it  was  decided  to  consider  the  effect  of  fluid 
capacitance  and  inertia.  The  connecting  line  was  modeled  by 
twc  resistance  blocks  separated  by  a  fluid  capacitance 
block.  The  breech  chamber  was  modeled  using  a  fluid 
capacitance  and  a  fluid  inertia  block,  as  was  done  for  the 
injector  chamber.  This  mcdel  gave  the  final  results 
described  in  the  section  comparing  computer  and  experimental 
data.  Ey  combining  these  building  blocks,  the  governing 
egtaticns  became:  (See  Fig.  (1)  for  notation) 


For  the  in*ectcr: 


1.  Pr  A  -PA  =  Np  V  '+  Kfs  V  +  Pdi  A 

111  1  11 


2.  V  -  V  =  ^P 

1  2  K  1 

3.  P  -  P  =  pL 

2  3/12 


For  the  connecting  line: 
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4.  P  -  P  =  Rv  V 

3  4  14 


5.  V  -  V  =  —  P 
4  5  K  4 


6.  p  -  p  =  Rv  V 

5  6  2  5 


For  th€  ireech  chamber: 


7. 


6 


8. 


K 


P 

8 


9. 


0 

P  A  -  Pds  A  =  Ms  V  +  Kfs  V 

8  3  ,  3  8  8 


These  ecuations  were  subsequently  modified  tc  account 
for  changes  in  cross-sectional  area  and  changing  geometry  as 
the  injector  chamber  decreased  in  volume  and  the  breech 
chamber  increased.  Once  the  projectile  reached  the  end  of 
the  breech  chamber,  the  force  balance  equation  was  nc  longer 
considered  applicable.  At  this  time,  it  was  noted  that  the 
remaining  eccations  for  the  connecting  line  and  the  breech 
chamber  cculd  be  combined  to  form  a  second  order 
differential  equation  for  the  breech  chamber  pressure  which 
described  the  "water  hammer"  type  pressure  oscillations 
displayed  by  the  experimental  data.  Thus, 
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Onc€  th€  injector  piston  reached  the  end  of  its  travel, 
the  injectcr  chamber  no  longer  existed  and  the  governing 
egnaticns  %ere  considered  as  no  longer  applicable. 
Dltimatelj,  the  injector  chamber  pressure  was  propagated 
down  the  connecting  line  as  the  propellant  fluid  came  to 
rest.  Throughout  these  analyses,  input  ram  pressure  was 
modeled  using  expcrential  terms  to  fit  experimental  data. 


The  eguaticns  frcm  the  lumped  parameter  approach  fcr  the 
propellant  fluid  were  combined  with  the  force  balance 
eguations  fcr  the  injector  piston  and  the  projectile  using 
State  Variable  methods.  Subsequently,  a  computer  program 
was  constructed  to  solve  these  state  variable  relations 
using  a  fourth  order  Runga-Kutta  integration  routine  for 
simultaneous  first  order  differential  eguations  which  was 
developed  at  the  Naval  Postgraduate  School.  A  brief 
description  cf  the  state  variable  method  and  a  listing  of 
the  ccmputer  program  can  be  found  in  Appendix  A. 

The  following  is  a  listing  cf  the  nomenclature  and  the 
values  cf  constants  used  in  the  derivation  cf  the  governing 
equations  anc  the  computer  program: 


A1 

Injectcr 

Cross-Sectional  Area 

1.77 

2 

in 

11. 42 

2 

cir 

A2 

Connecting  Line 

Cross-Sectional  Area 

.255 

2 

in 

1.65 

2 

car 

A3 

Ereech  Chamber 

.49 

2 

in 

3.16 

2 

cs 

Cross-Sectional  Area 

iR 

Ram  listen 

Cross-Sectional  Area 

1.77 

2 

in 

11.42 

2 

cm 

11 

Injectcr  Length 

1.60 

in 

4.22 

cir 
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L2 

Connecting  Line  Length 

16.0 

in 

40.6 

cm 

L3 

Breech  Chamber  Length 

5.0 

in 

12.7 

cm 

MF 

Injectcr  Piston  Mass 

.0052 

2 

lb~sec 

913.8 

gm 

in 

MS 

Projectile  Mass 

.00053 

2 

lb~sec 

93.  1 

gm 

in 

KFP 

Injectcr  Frictipn 

.001 

lb-sec 

175.7 

gm/sec 

Factor 

in 

KFS 

Projectile  Friction 

.001 

lb-sec 

175.7 

gm/sec 

Factcr 

in 

FDI 

Injector  Back  Pressure 

24.0 

PSI 

163.3 

kFa 

FES 

Projectile  Back 

.50 

PSI 

3.4 

kFa 

Pressure 

K 

Effective  Bulk  Modulus 

3200 

PSI 

6808 

kFa 

fiHC 

Propellant  Density 

.000093 

2 

lb-sec. 

.994 

3 

gu/cm 

4 

in 

FB 

fiam  Pressure 

Computed 

FI 

Injectcr  Chamber  Pressure 

Computed 

VI 

Injectcr  Piston  Velocity 

Computed 

V2 

Injector  Chamber  Exit  Velocity 

Computed 

P4 

Connecting  Line  Pressure 

Computed 

V6 

Breech  Chamber  Entrance 

Velocity 

Computed 

F8 

Breech  Chamber  Pressure 

Computed 

V8 

Projectile  Velocity 

Computed 

Th€  treatttent  cf  several  areas  of  the  computer  model 
continually  reappeared  as  requiring  refinement.  The 
governing  differential  equations  worked  well  during  the 
dynamic  pcrticn  cf  the  feed  cycle  but  experienced 
difficulties  during  the  initial  and  final  static  periods. 
The  handling  cf  static  friction  and  back  pressure  as  a 
ccrstant  value  created  the  possibility  of  negative 
velocities  urtil  the  driving  ram  pressure  overcame  the 
system  tack  pressure.  These  negative  values  never  occurred 
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in  the  real  system  due  to  the  geometric  restraints  on  the 
injector  piston  and  projectile;  therefore,  the  computer 
model  hac  tc  be  manipulated  to  maintain  this  condition. 
Unfortunately  the  need  for  simultaneous  solution  cf  the 
governing  eguations  made  it  difficult  to  manipulate  the 
initial  ccnditions  without  greatly  increasing  the  complexity 
of  the  computer  modeJL.  After  the  projectile  stopped,  the 
transitici  from  the  force  balance  equation  on  the  projectile 
to  an  equation  describing  the  fluid  velocity  in  response  to 
the  pressure  transients  was  awkward.  No  simple  differential 
equation  describes  the  complexity  of  the  wave  mechanics 
involved  with  the  reflection  of  the  pressure  waves  in  the 
system.  In  this  regard  several  alternate  approaches 
involving  sequential  alterations  of  the  governing  eguations 
were  tried  with  varying  degrees  of  success. 
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description  CF  experimental  APFAMIi^S 


In  cid€i  that  the  comparison  of  the  experimental  and 
ccuputer  generated  data  can  be  fully  understood,  a  brief 
description  cf  the  NPS  experimental  apparatus  and  the 
conduct  cf  tie  associated  experiments  is  included  (See  Ref. 
1 ,  pages  17-27) . 

The  test  chamber  was  fabricated  frcm  a  three-inch  O.D. 
Lucite  cjlirder,  18  inches  long,  bored  to  a  20  mm  inside 
diameter  and  fitted  with  aluminum  end  caps.  The  chamber  was 
loaded  with  e  brass  slug  weighing  93  grams  which  rode  cn  two 
graphite  filled  Teflpn  sealing  rings.  The  brass  slug,  which 
simulated  the  projectile,  was  cycled  from  the  breech  end  of 
the  chamber  tc  the  barrel  end  and  returned  to  the  breech 
end,  completing  one  hypothetical  firing  cycle. 

Because  cf  the  desire  to  vary  the  charge  to  mass  ratio, 
a  variable  chamber  velocity  was  necessary.  To  accomplish 
this  with  one  chamber,  a  volume  control  retaining  rod  was 
designed  into  the  system.  This  brass  rod,  bored  to  allow 
gas  to  pass  its  length,  was  threaded  through  a  plate  which 
was  attached  to  the  barrel  end  cap  holding  the  rod  in  the 
chamber.  The  rod,  which  has  a  Teflon  disc  cn  the  end,  not 
cnly  established  the  volume  of  the  test  chamber,  but 
provided  a  buffer  stop  for  the  slug  at  the  end  of  its 
forward  motion.  Another  Teflon  buffer  was  affixed  to  the 
breech  eud  cap  to  cushion  the  slug  in  return  motion. 

The  LEG  simulator  started  a  simulated  firing  cycle  with 
the  slug  at  the  breech  end  of  the  empty  chamber,  as  shewn  in 
Fig.  (3). 
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Power  Piston  Injector  P  “i  Test  Chamber 
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Figure  3  -  PNEUMATIC  CONTROL  SYSTEM  IN  RAM  POSITION 


The  simulated 


This  is  the  ready-to-ram  position, 

propellant,  distilled  water,  was  then  introduced,  ramniing 
the  slug  to  the  opposite  (barrel)  end  of  the  chamber  as  the 
chamber  was  filled.  This  was  accomplished  by  applying  gas 
pressure  tc  a  power  piston  which  drove  the  injector  piston. 
The  injector  piston  forced  the  propellant  past  a  flow  check 
valve  and  into  the  chamber.  This  placed  the  slug  in  the 

ready-tc-f ire  cr  ic-battery  position.  In  an  actual  gun,  the 
prcpellant  wculd  be  ignited  at  this  time  in  the  cycle.  Due 
to  laboratory  constraints,  an  expulsion  system  was  used. 

An  IVDT  displacement  transducer  was  manufactured  and 
mounted  next  tc  the  injector  piston.  The  LVDT  was  attached 
tc  the  connecting  rod,  between  the  power  piston  and  the 
injector  pistcn.  The  volume  of  the  liquid  being  placed 
under  pressure  during  each  shot  was  measured  by  filling  the 
system  in  the  ready-to-ram  position  and  then  draining  it 
into  a  graduated  beaker.  By  measuring  the  displacement  of 
the  injector  piston  head  with  the  LVDT,  the  volumetric  rate 
of  fluid  injection  during  the  ram  stroke  was  obtained. 

Two  pressure  taps  were  drilled  in  the  test  chamber. 
These  taps  were  located  as  close  to  the  breech  end  as 
possible,  cne  at  2Cr  degrees  frcm  top  center  and  the  other  at 
20  degrees  from  bottom  center  in  a  counterclockwise 
direction,  as  seen  from  the  breech  end  of  the  chamber.  Only 
the  bottCE  Iccation  was  used  to  record  data. 

A  4-chanrel  Hewlett-Packard  3960  Magnetic  Tape  Recorder 
was  used  to  IM  record  desired  data  during  system  operation. 
For  each  cata  run  signals  were  recorded  from  two  Kaman 
diaphragms  type  ;(1000  psi)  pressure  transducers,  one 
connected  to  the  breech  pressure  tap  on  the  test  chamber, 
and  the  ether  to  the  gas  injection  side  of  the  power  piston. 
The  pressure  signals  were  processed  with  a  Kaman  Digi-Vit 
Readout  Dnit  which  also  provided  a  visual  (digital)  display. 
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Data  froff  the  IVDT  displacement  transducer  was  also  recorded 
cn  the  tape, 

A  Brush  Becorder  (Mark  280)  was  used  to  obtain  a  visual 
display  cf  the  recorded  data.  By  transcribing  the  desired 
signals  cn  the  Magnetic  Tape  fiecorder  at  a  tape  speed  cf  15 
feet  per  minute  and  playing  them  back  into  the  Brush 
Becorder  at  3-3/4  feet  per  minute,  the  time  scale  of  the 
output  was  expanded  by  a  factor  of  four  on  the  Brush 
recordings  (viz.,  from  a  real-time  maximum  cf  200  mm/sec  to 
a  delayed  tine  maximum  of  800  mm/sec.) 
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V 


CCKPARIS^  OF  COMPOTEB  AND  EXPERIMENTAL  DAjlA 


As  iD€iiticEed  previously,  experimental  data  was  taken  for 
ran  gas  pressures  ranging  from  50  to  220  psig,  Fron  this 
ccllecticn  of  data,  pne  run  at  140  psig  was  selected  as 
representative  of  system  performance.  It  was  felt  that  any 
model  which  vculd  adequately  describe  system  operation  at 
this  intermediate  pressure  would  be  valid  for  the  entire 
range  of  expected  IPG  driving  pressures. 

Figure  (^l)  is  a  comparison  cf  computed  and  recorded  data 
for  injectcr  pistor  displacement.  As  can  be  seen,  the 
agreement  is  very  good.  This  is  not  very  significant  in 
that  all  models  tried,  as  well  as  the  original  analog 
computer  model,  were  able  to  correctly  predict  injector 
pistcn  metier. 

Figure  (5)  is  a  comparison  of  analytical  results  and 
experimental  data  for  breech  chamber  pressure.  The  computer 
model  fellows  the  shape  of  the  experimental  curve  fer  the 
duration  cf  the  time  that  the  projectile  is  in  motion  (0-30 
msec) .  It  oscillates  very  rapidly  but  does  not  fall  eff  to 
a  zero  value  at  the  end  of  the  projectile  motion.  The 
frequency  of  these  initial  oscillations  is  approximate! j  1.0 
khz  which  is  within  the  -frequency  range  of  the  Kaman 
pressure  transducers  used  in  testing.  However,  the  mounting 
of  these  transducers  within  a  connecting  cavity  instead  of 
flush  with  the  chamber  wall  could  have  led  to  them 
experiencing  a  reduced,  lagging  frequency  response.  It  is 
felt  that  seme,  if  not  all,  of  the  computed  pressure 
oscillaticns  must  exist  as  evidenced  by  the  close 
correspondence  of  the  first  two  peaks  in  Fig.  (5) .  The 
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accuracy  cf  the  experimental  pressure  reading  could  also  be 
coDsidered  as  being  reduced  by  the  location  of  the  pressure 
transducer  in  the  test  setup.  It  is  possible  that  the 
monitcring  fiessure  tap  was  too  close  to  the  end  cf  the 
breech  chamber  and  may  not  have  sensed  the  full  chamber 
pressure  during  the  latter  part  of  the  projectile’s  travel 
when  the  fluid  velocity  is  greatest.  It  is  felt  that  to 
adeguately  describe  the  pressure  decrease  toward  zero,  it 
would  be  necessary  to  completely  account  for  exact  changes 
in  system  cecmetry.  To  do  this  would  require  using  a 
distributed  parameter  approach  featuring  variable  lengths 
fcr  the  injector  chamber  and  breech  chamber. 

The  pressure  transients  which  occur  after  the  projectile 
steps  are  ccnplex  interference  phenomena  which  are  not  fully 
described  by  the  mcdel.  However,  the  peak  pressure  and  the 
natural  freguency  of  oscillation  of  the  computer  model 
cempare  fevcrably  with  the  experimental  data. 
Onf ortunatel j ,  the  damping  characteristics  cf  the  model  do 
not  follow  the  experimental  data  well.  The  system 
parameters  which  determine  these  values  will  be  discussed  in 
the  next  section. 

It  should  be  '  reiterated  that  the  lumped  parameter 
approach  is  at  approximation..  In  a  distributed  parameter 
approach,  pressure  and  velocity  would  vary  continuously  with 
distance  within  the  system  as  well  as  with  time.  In  the 
lumped  paraneter  approach,  these  changes  are  assumed  to 
occur  only  at  the  input  and  output  of  a  building  block; 
therefore,  pressure  and  velocity  are  considered  constant 
within  the  building  block.  This  will  always  lead  to  some 
discrepancies  when  comparing  model  data  to  experimental  data 
taken  at  at  a  fixed  point. 

Figure  (6)  shows  the  simultaneous  pressure  history  for 
several  areas  of  the  LPG  feed  system.  The  input  ram 
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pressure,  injectcr  chamber  pressure,  connecting  line 
pressure,  ard  breech  chamber  pressure  are  graphed  as 
furcticns  cf  time. 

Figures  {7)  and  ,(8)  show  the  simultaneous  velocities  at 
several  pcirts  as  predicted  by  the  computer  model.  The 
irjectcr  pistcn  velocity,  fluid  propellant  velocities  at  the 
beginning  and  end  of  the  connecting  line,  and  projectile 
velocity  are  plotted  as  functions  of  time. 

The  injector  and  projectile  displacements  as  functions 
of  time  are  included  as  Figure  (9) .  As  can  be  seen  from 
Figures  (7) ,  (8) ,  and  (9) ,  the  movement  of  the  injector 
piston  is  much  mere  rapid  than  that  of  the  projectile.  It 
is  felt  that  this  velocity  difference  is  very  much  a 
function  cf  system  geometry  and  should  not  be  considered  as 
a  generali2ation  for  all  LPG  feed  systems.  The  NPS 
experimental  apparatus  being  modeled  had  only  small  pressure 
drops  between  injector  and  breech  chambers.  In  addition, 
tfce  projectile  slug  was  an  order  of  magnitude  lighter  than 
the  injector  piston.  In  large  scale  systems,  with 
significant  pressure  drops  and  very  massive  projectiles,  it 
is  quite  possible  that  the  injector  piston's  full  stroke 
will  occur  significantly  before  the  projectile  has  seated. 
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Figure  4  -  COMPARISON  OF  INJECTOR  PISTON  DISPLACEMENT 
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Figure  5  -  COMPARISON  OF  BREECH  CHAMBER  PRESSURE 
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Time  (Sec.  x  10^) 

Figure  6  -  SIMULTANEOUS  PRESSURE  HISTORIES 
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Figure  7  -  SIMULTANEOUS  INJECTOR 

VELOCITY  HISTORIES 
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Figure  8  -  SIMULTANEOUS  BREECH  CHAMBER  VELOCITY  HISTORIES 
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VI 


MALISLS  OF  SXSTEM  PERFC5MANCE 


A.  PABAfiETEfS  AFFECTING  PHOJECTILE  RAN  TIME 

Several  parameters  can  be  varied  to  "tune”  the  model  to 
correct  ram  time.  In  equations  (1)  and  (9)  (the  force 
balances  cn  the  injector  piston  and  projectile)  empirical 
constants  are  included  which  account  for  any  static  friction 
and  back  pressure  (PDI  =  injector  back  pressure,  PDS  = 
projectile  lack  pressure).  In  equations  (4)  and  (6), 
resistance  ccefficients  are  introduced  to  account  for  system 
resistances.  Finally  in  equations  (2),  (5),  and  (8), 
propellant  effective  bulk  modulus  values  appear. 

The  value  for  the  injector  and  projectile  total  tack 
pressure  is  difficult  to  determine  precisely.  Reference  1 
cited  an  experimentally  determined  value  of  24  psi.  The 
analog  computer  model  presented  in  Ref.  2  used  a  value  of  46 
psi  for  a  ran  pressure  of  140  psi.  The  present  model  uses  a 
value  of  24  psi  for  PDI  and  0.5  psi  for  PDS,  or  a  total  of 
24.5  psi  tack  pressure,  which  agrees  with  the  experimental 
value.  The  value  of  system  back  pressure  will  vary  with 
driving  pressure  and  the  configuration  of  each  LPG  system's 
injector  and  kreech  chambers. 

Some  basic  research  has  been  conducted  to  relate  the 
value  of  the  resistance  ccefficients  to  fluid  properties. 
Unfortunately  these  studies  which  are  described  in  Refs.  4 
and  5,  dealt  with  steady  state  fluid  flow  in  constant 
diameter  lices  without  flow  restrictions ,  making  their 
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results  inapplicable  to  the  present  LPG  feed  system  prcblem. 
The  value  cf  BV^  can  be  related  to  the  damping  ratio  of  the 

"water  hammer"  pressure  oscillations.  The  value  of  was 

determined  ty  an  iterative  process  to  obtain  the  best  fit  to 


experimental  data.  Figure 

ram-time  on  BV  and  RV  . 

1  2 


(10)  shows  the  dependence  of 
keeping  all  ether  variables 


constant.  is  expected  intuitively,  increasing  the  fluid 


resistance  creates  a  larger  system  pressure  drop,  resulting 
in  a  lower  pressure  exerted  on  the  face  cf  the  projectile 
and  hence  slower  ram  times.  The  slope  cf  the  two  curves  are 
almost  identical  sp  that  no  significant  advantage  would  be 
achieved  by  system  designs  which  try  to  minimize  either  RV 


or  RV  to  the  detriment  of  the  other. 
2 


It  should  be  noted. 


however,  that  higher  values  for  RV^  do  tend  to  slow  down  the 
ram  time  acre  than  high  values  for  RV  .  Since  the  magnitude 


of  BV^  would  probably  depend  mostly  on  the  pressure  drep  at 
the  gun  valve,  which  seals  the  breech  chamber,  its  design 


should  he  closely  watched  to  ensure  rapid  ram  times. 
Likewise  the  system  designer  will  have  to  pay  close 
attention  to  the  design  of  the  piping  and  valves  in  a  large 
scale  LPG  feed  system  to  achieve  optimum  performance. 


The  value  cf  the  effective  bulk  modulus  was  determined 
from  the  undamped  natural  frequency  of  the  "water  hammer" 
pressure  cscillaticns  as  will  be  discussed  in  the  next 
section.  The  bulk  modulus  is  a  fluid  parameter  which 
characterizes  the  spring  effect  of  a  liquid.  The  bulk 
modulus  can  be  substantially  lowered  by  the  elasticity  of 
the  chamber  and  connecting  line  walls  and  the  ameunt  of 
entrapped  gas  present  in  the  propellant  fluid.  As  can  be 
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S€€D  froff  Figure  (11) ,  the  value  of  the  bulk  modulus  has 
coly  a  small  effect  on  ram  time.  Between  bulk  modulus 
values  of  10,000  and  100,000  psi  (the  expected  operating 
region  of  an  operational  system)  the  ram  time  is  almost 
ccEStant,  varying  less  than  a  guarter  of  a  millisecond. 
This  is  an  eccouragijig  result  since  minimizing  entrapped  gas 
is  a  difficult  and  cpstly  design  constraint. 
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Figure  10  -  RESISTANCE  COEFFICIENTS  VS  BAM  TIME 
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E.  PAEAMiTEES  AEEECTING  CHAMBER  PRESSORE  CSCILLATICNS 

As  eiplained  in  Section  III,  after  the  projectile  has 
stepped  the  chamber  pressure  can  be  described  by  a  second 
order  differential  equation  (Egn.  10) .  If  the  system 
described  by  this  differential  equation  is  underdamped,  the 
solution  to  the  equation  is  an  exponentially  damped  sinusoid 
with  a  characteristic  damping  ratio  and  undamped  natural 
frequency. 

From  Egn.  (10),  the  system  damping  ratio,  ^  is: 


A 

A2  2y/?  K 


If  this  daiiping  ratio  is  increased,  the  peak  chamber 
pressure  as  well  as  the  dissipation  period  of  the  system 
pressure  oscillations  will  be  decreased.  From  a  designer's 
point  of  ^iew,  it  would  be  advantageous  for  both  of  these 
quantities  to  be  as  low  as  possible.  Hence,  it  would  seem 
that  an  increase  in  system  resistance  to  enhance  damping 
would  be  desirable.  However,  as  shown  in  the  previous 
section, an  increase  in  the  system  resistance  also  increases 
the  ram  time-  an  undesirable  occurrence.  Thus,  alterations 
in  system  configuration  which  affect  resistance  coefficient 
values  will  have  to  be  accomplished  carefully  tc  insure 
optimal  sjstem  performance. 

From  Egn.  (1C) #  the  undamped  natural  frequency  cf  the 
pressure  cscillations  in  the  LPG  system  is: 


The  primary  effect  of  ullage  would  be  to  decrease  the  value 
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of  effective  hulk  modulus  and  conseguently  tc  decrease  the 
natural  tregoency  of  pressure  oscillations.  Hence,  should 
it  be  found  in  future  tests  that  the  actual  liquid 
propellant  is  highly  sensitive  to  pressure  oscillations, 
changes  in  ullage  may  be  effected  to  decrease  system 
ringing.  Fortunately,  as  shown  in  the  previous  section,  the 
effect  of  ullage  variations  on  ram  time  is  minimal. 
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VII.  CONCLGSIONS 


Several  important  conclusions  can  be  drawn  frcm  the 
output  of  the  IPG  ram  feed  computer  model.  Thus^  it  is 
realized  that  serious  design  consideration  must  be  given  to 
the  values  of  fluid  resistances  in  the  feed  system  in  order 
tc  minimize  ram  time  and  peak  pressures^  and  to  optimize  the 
damping  cf  system  pressure  transients.  Fortuantely^ 
ccncommitant  model  results  also  show  that  the  effect  of 
entrapped  gas  is  minimal  for  the  overall  operation  of  a  LPG 
feed  system.  The  model  was  unable  to  account  fcr  any 
cavitation  effects  so  this  is  an  area  that  must  be  a  subject 
cf  further  studies. 

The  fact  that  large  pressure  transients  cccur  during  the 
projectile  ramming  cycle  indicates  the  difficulties  that 
might  arise  with  proposed  designs  for  liquid  propellant  guns 
with  variable  propellant  volumes  and  projectile 
displacements .  The  pressure  transients  experienced  in  these 
systems  will  undoubtedly  be  very  complex  in  nature. 

If  further  improvements  are  desired  in  the  modeling  of 
LPG  feed  systems^  the  next  step  should  be  the  inclusion  of 
thermody ramie  effects.  For  any  real  fluid,  resultant  feed 
system  pressure  and  velocity  changes  will  cause  temperature 
changes  which  may  substantially  affect  the  propellant  fluid 
density.  Consequently,  propellant  ignition  characteristics 
at  high  rates  cf  fire  will  undoubtedly  be  affected. 

The  lumped  parameter  LPG  feed  system  computer  model  has 
been  shown  tc  be  in  favorable  agreement  with  experimental 
data.  It  is  felt  that  this  model  has  been  sufficiently 
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validiated  tc  allow  its  use  in  more  complicated  LEG  feed 
system  designs.  It  is  expected  that  this  model  will  be  used 
by  the  Naval  Ordnance  Station,  Indian  Head,  Maryland  to 
assist  in  studying  the  fluid  dynamic  characteristics  of  a  30 
mm  scale  model  LPG  feed  system  becoming  operational  in  June 
1S76.  It  will  be  necessary  to  make  the  obvious  changes  to 
describe  the  different  geometry  of  the  30  mm  scale  model. 
The  additional  line  lengths,  and  the  larger  pressure  drops 
due  to  complex  valve  arrangements  in  the  30  mm  system  will 
require  different  values  for  model  back  pressure  and 
resistance  coefficients.  The  presence  of  an  accumulator 
near  the  breech  chamber  will  affect  the  value  of  the  system 
damping  ratio.  It  is  possible  that  the  different  geometry 
and  mass  characteristics  of  the  30  mm  model  could  result  in 
the  injector  piston  being  slower  than  the  projectile. 
However  by  incorporating  all  of  these  alterations,  the 
manner  in  which  each  of  the  system  parameters  affect  total 
system  performance  will  be  able  to  be  predicted  by  the 
versatile  lumped  parameter  computer  model  developed  in  this 
study. 

As  such,  computer  simulation  is  a  useful  design  tool. 
However,  tc  remain  useful,  it  must  be  viewed  in  its  proper 
perspective.  Full  scale  prototype  testing  is  the  only 
conclusive  method  of  demonstrating  ordnance  system 
performance.  Unfortunately,  testing  is  both  costly  and  time 
consuming.  Thus,  computer  simulation,  no  matter  how  simple 
or  complex,  can  be  used  to  designate  critical  testing 
instances  and  to  identify  those  areas  where  design  efforts 
will  be  icst  productive.  In  this  way,  the  time  involved  in 
optimizing  system  performance  as  well  as  subsequent 
prcducticn  costs  can  be  reduced.  In  explosives  research  and 
ordnance  design,  which  is  still  an  empirical  and  necessarily 
hazardous  science,  computer  simulation  can  be  particularly 
useful  in  specifying  the  correct  approach  for  design  testing 
practices.  In  this  respect,  such  models  as  the  lumped 
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parametei  model  developed  in  this  study  can  fce  considered  as 
critical  signposts  at  some  of  the  many  crossroads  in  systems 
development.  However,  it  must  be  realized  that  as  such  they 
cnly  point  the  way.  The  vehicle  for  arriving  at  an 
operational  system  can  only  be  diligent  research  and 
engineering  based  on  a  measured  progression  of  system 
demcnstraticns  overseen  by  dedicated  project  managers  aware 
of  the  many  pitfalls  along  the  way. 
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APPENDIX  A 


COMPOTER  PROGRAM  LISTING 


The  hasis  of  the  state  variable  method  of  systems 
analysis  is  the  interpretation  that  the  energy  state  of  the 
cciponents  of  a  system  completely  describes  the  condition  of 
a  system.  As  defined  in  Ref.  6,  the  state  of  a  system  is  the 
set  of  variables,  the  "state  variables",  which  contain 
sufficient  information  about  the  present  condition  of  the 
system  tc  permit  the  determination  of  all  future  time 
history  cf  the  system  -  provided  that  all  future  inputs  to 
the  system  are  known.  Therefore  the  energy  state  cf  those 
elements  which  store  energy  completely  describes  the  system. 
In  the  case  of  the  lumped  parameter  model,  these  elements 
are  the  fluid  inertia  and  capacitance  and  the  inertia  cf  the 
injectcr  pistcn  and  projectile.  The  energy  state  cf  these 
storage  elements  can  be  described  as  a  function  of  time  in 
terms  of  the  state  variables  ->•  pressure,  velocity,  and  the 
system  input  pressure.  Only  those  variables  required  to 
completely  specify  the  state  of  the  system  need  to  be 
included • 

The  first  step  in  arriving  at  the  computer  program  was 
to  take  equations  (1>  through  (9)  and  normalize,  or  solve 
them  for  the  highest  derivative.  The  state  variables  then 
become  the  pressure  or  velocity  associated  with  these  first 
derivatives.  The  computer  program  was  developed  by 
ccmbining  equations  and  defining  new  variables  as  follows: 

X(1)  = 
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X(2)  =  E 
X(3)  =  V. 
x(4)  =  e;; 

X(5)  = 

X(8)  =  Eg 
X(9)  =  Eg 

Eguatioo  (1)  becones: 

-Kfp  At  Ar  At  „ 

XDOTQ) - XQ)  -  --  X(2)  +  --  Pr  -  --  Pdi 

Mp  Mp  Mp  Mp 

In  the  ccmputer  proigram,  PR  is  represented  ty  X(30)  and  KFP 
and  PDI  heccme  program  constants,  C <7)  and  C  (2) , 
respectively.  The  value  of  C(2)  and  C(7)  are  input  on  data 


cards. 


Osing  the  fact  that  P^=P^,  equation  (2)  becomes: 


•  XD0T(2)  =  X(l)  -  -  X(3) 

Li 

where  K  is  input  to  the  computer  program  as  C(3).  By 

ccmbining  equations  ,(3)  and  (4)  and  using  the  incompressible 

flow  assuapticn  A  V  =A  V  the  next  equation  is  derived: 

12  2  4 


XD0T(3)  = 


Rho  L 


(  X(2)  -  Rv  ^  X(3)  -  X(4)  ) 


where  HV  becomes  C  (4) .  Next 
1 

incompressible  flow  assumptions 
yield: 


equation 


A  V  =A  V 
12  2  4 


(5) 

and 


and 


the 


A  V  =A  V 
2  5  3  6 


At  K  A,  K 

XD0T(4)  =  —  —  X(3)  -  -  -  X(4) 

^2  L2  A2  L2 

and  as  before,  C(3)=K.  Equations  (6)  and  (7)  combine  with 

the  flow  velocity  equation  A  V  =A  V  to  yield: 

2  5  3  6 


XD0T(5)  = 


A3 


where  E7 


_  (  X(«  -  ^  X(5)  -  X(8)  ) 

becomes  C  (5) .  Since  V  =V  ,  equation  (8)  becomes: 

6  7 


XD0T(8)  =  “  X(5)  -  p-  X(9) 
La  La 
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where  K  again  beccmes  C(3),  Finally,  equation  (9)  fceccnies: 

■“Kfs  a  Ao 

XD0T(9)  =  -  X(9)  +  X(8)  -  Pds 

Ms  Ms  .  Ms 

where  KFS  and  FDS  become  C  (8)  and  C  (9)  respectively.  The 
auxiliary  equations: 

XD0T(6)  =  X(l)  ■ 

XD0T(7)  =  X(9) 

are  used  tc  calculate  injector  piston  displacement,  X  (6) , 
and  projectile  displacemen't,  X(7). 
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IT  IS  ALSO  POSSIBLE  TO  PROCESS  THE  INTEGRATION  IN  EITHER  TVO 
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